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PATENT 



METHOD FOR INCREASEVG ACCURACY OF 
MEASUREMENT OF MEAN POLARIZATION MODE DISPERSION 

Field of the Invention 

[0001] The invention relates to the field of high-speed optical telecommunication systems, 
and more particularly to a method for increasing the accuracy of the measurement of mean 
polarization mode dispersion in optical fiber. 

Background of the Invention 

[0002] The optical waveguides or fibers used to transmit signals in optical 
telecommunication systems are characterized, in part, by the vector property of polarization 
mode dispersion (PMD). Polarization mode dispersion occurs as a result of birefiingence in 
the fiber, which may be caused by physical asymmetry in the fiber construction itself, or by 
stress, stram, or other extemal forces imposed on the fiber. In addition, random polarization 
coupling can occur, giving rise to a time- varying statistical factor. Optical fibers display an 
anisotropy in the refiractive index, which v^ll vary as a function of position and time. 
Consequently, components of an optical signal that differ in polarization will propagate at 
different velocities, resulting in a differential group delay (DGD) between the components, 
and causing significant broadening of the optical pulses propagating along long lengths of 
fiber. 



[0003] The PMD is fully characterized by a vector quantity r(fi?) where the DGD is the 

magnitude of the vector t(o}) . As shown in FIG. 1, the DGD is generally designated as r 

10. Any state of polarization (SOP) can be resolved into directional components along two 
orthogonal principal states of polarization (PSP) 12, 14. The DGD or r 10 then represents 
the separation in time between fast PSP 12 and slow PSP 14, after traversing a length of 
optical fiber 16. For each optical firequency or wavelength propagating in a fiber, there 
always exists two PSPs, such that the pulse spreading due to the first-order PMD vanishes if 
only one PSP is excited. The PMD is typically characterized in terms of an average DGD 
corresponding to different frequencies, and is independent, to first order, of wavelength, 
temperature, and extemal perturbations. In low mode coupled fiber, this measure of DGD 



averaged over a large range of optical frequencies is fairly constant over time, but in high 
mode coupled fibers, for example, in long fiber spans, the fi-equency-averaged DGD varies 
randomly in time, due to the combined effects of the variations in birefiingence and random 
polarization mode coupling along the fiber length. This statistical variation in DGD lends 
itself to characterization of the DGD in terms of a statistical figure of merit, mean DGD. 

[0004] Higher orders of the polarization mode dispersion also exhibit statistical properties. 
The effect of second-order polarization mode dispersion (SOPMD) 18 is shown in FIG. 1. 
The SOPMD is the first derivative of the PMD with respect to fi-equency, representing the 
change in the PMD as a fimction of fi-equency. The SOPMD, therefore, additionally 
characterizes the overall pulse spreading due to the firequency-dependence of the PMD and 
the spectral bandwidth of the injected optical pulse 19. 

[0005] The polarization mode dispersion of a fiber is imlike most other sources of 
degradation in an optical telecommunication system, in its dependence on both time and 
fi*equency. Conventional methods for characterizing the fiiU PMD vector over a fi-equency 
range, well known by those skilled in the art, include the Poincare Sphere Analysis (PSA), 
the Jones Matrix Eigenvalue (JME), MuUer Matrix Method (MMM), Fixed Analyzer and 
interferometric techniques. These methods provide a measure of mean DGD and root mean 
square (RMS) DGD, which is calculated from the set of frequency-dependent DGD values. It 
is then commonly assumed by those skilled in the art that the statistical DGD follows a 
Maxwellian distribution, so that a true mean DGD (r) , determined by averaging the DGD 
values obtained for a nvimber of fibers over a bandwidth B of optical frequencies, can be 
estimated by multiplying the measured RMS DGD ^J{^^)~ by a factor of 

[0006] The fimdamental problem in accurately evaluating a statistical limitation to an 
estimation of the mean DGD of a fiber, in order to find a more precise measurement of the 
mean DGD, was first recognized in a paper by N. Gisin, B. Gisin, J.P. Von der Weid, and R. 
Passy, entitled "How Accurately Can One Measure a Statistical Quantity Like Polarization- 
Mode Dispersion?" IEEE Photon. Tech. Lett., Vol. 12, pp. 1671-1673 (August 1996), which 
is incorporated herein by reference. The accuracy of mean DGD estimation does improve as 
the mean is taken over a larger spectral bandwidth (approaching the ideal theoretical case 
where 5-^oo). However, contrary to the statistical requirement that each of the measurements 
used to calculate an average be independent, the DGD at nearby wavelengths are not 

-2- 




frequency independent. Gisin et al. demonstrated that this frequency dependence resulted in 
lower uncertainty in the mean DGD (around 9%) for larger PMD on the order of 1 
picosecond (ps) e.g., as compared to a 28% uncertainty in mean DGD measurement when the 
PMD is smaller (on the order of . 1 ps). The uncertainty in mean DGD measurement 
increases with decreasing source bandwidth. Gisin et al. demonstrated that the same level of 
uncertainty is intrinsic to all measurement techniques that average the DGD over wavelength. 

[0007] The mathematical formalism was developed further by M. Shtaif and A. Mecozzi, 
"Study of the Frequency Autocorrelation of the Differential Group Delay in Fibers with 
Polarization Mode Dispersion," IEEE Photon, Tech, Lett, Vol. 25, pp. 707-709 (May 2000), 
which is incorporated herein by reference. In measurements of the frequency autocorrelation 
of the DGD, the square DGD, and orientation of the PMD vector, Shtaif et al. showed that all 
corresponding correlation bandwidths are comparable. Shtaif et al. also showed that all 
statistical properties of the PMD characterizing the fiber imder test are uniquely defined by 
the mean DGD. 

[0008] Polarization mode dispersion (PMD) is recognized as a potentially limiting 
impairment for high-speed long-haul optical transmission. Moreover, precise measurement 
of the true mean differential group delay (DGD) of individual fiber links and whole fiber 
routes is important for accurate estimation of service outage probabilities. Since PMD varies 
with time, as well as with frequency, measurements of the mean frequency-averaged DGD of 
the same fiber taken at different times may differ from each other and from the true value of 
mean DGD for a given fiber. For DGD values in the usual range of interest, and within the 
optical bandwidths of commercially available equipment, the variance of DGD measurements 
is approximately inversely proportional to the optical bandwidth of the optical source used 
for the measurement. In other words, an accurate measurement of the mean DGD of low 
birefringence fiber is limited by the optical bandwidth of the source used for the 
measurement. 

[0009] The need for precise PMD characterization will increase as the high-speed networks 
of the fiiture employ very low PMD fibers. There exists a need, therefore, for more precise 
measurement of the mean DGD of individual fiber links and whole fiber routes than is 
presently provided by conventional methods. 
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Summary of the Invention 

[0010] The present invention, which addresses the needs unmet by conventional methods, 
relates to methods of improving accuracy of measuring a differential group delay (DGD) in 
an optical fiber link and in a whole optical fiber route. 

[001 1] A method of the present invention for measuring a tme mean differential group delay 
(r) of at least a length of optical fiber includes the initial step of measuring a mean square 

differential group delay (^^)^ averaged over a finite spectral bandwidth 5 of a source, using 
a polarization mode dispersion measurement apparatus. A root mean square differential 
group delay is then calculated in accordance with ^^r^^^ and a systematic correction factor 
8 to the conventional method of estimating true mean {r) fi-om the measured root mean 
square differential group delay ^(r^)^ is applied. The systematic correction factor e 
minimizes a systematic error caused by the finite spectral bandwidth of the source. 

[0012] Preferably, the systematic correction factor e is applied to the mean square differential 
group delay in accordance with (r) = ^J_^^2^^ + to obtain the true mean differential 

group delay, (r). In the regime where r^^B » 1, (r) is calculated according to: 

(16a) and (16b) 

8 1 



in other words, e is 



9^/2 B 



[0013] This method, which is applied directly to a measured mean square differential group 
delay, may be applied to measurements taken using time-domain techniques with an 
apparatus such as an interferometer. This method may also be applied to measurements taken 
using frequency-domain techniques, such as Jones Matrix Eigenanalysis, Poincare Sphere 
Analysis, and MuUer Matrix Method using an apparatus, for ex£imple, including a 
polarimeter. The method may be used for measuring the mean differential group delay 
through a single optical fiber link, or an entire optical fiber route. 
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[0014] In another embodiment of the method of the present invention, a method for 
measuring a mean differential group delay (r) of at least one length of optical fiber, includes 

an initial step of characterizing a polarization mode dispersion vector as a function of 
frequency using a frequency-domain polarization mode dispersion measurement apparatus. 
The method further includes calculating a second-order polarization mode dispersion vector 
as a function of frequency from the polarization mode dispersion vector, and calculating a 
mean of the square root of a magnitude of the second-order polarization mode dispersion 

vector |f ^ I to obtain a first result, according to l^f^ | ^ ^ . The first result is multiplied by a 

proportionality coefficient ^2 to calculate the mean differential group delay, in accordance 
with the following equation: 



(KI^) = (->- (21) 



[0015] Preferably, the proportionality coefficient ^2 is obtained from the probability density 
function of the second-order polarization mode dispersion vector. Most preferably, ^^2 is 
substantially equal to 1.37. 

[0016] A further embodiment of the method of the present invention provides a method for 
measuring a mean differential group delay (r) of at least one length of optical fiber, 

including an initial step of measuring a magnitude of a polarization mode dispersion vector as 
a function of frequency, using a frequency-domain polarization mode dispersion 
measurement apparatus, where the magnitude of the polarization mode dispersion vector is a 
scalar differential group delay. The method further includes calculating a frequency- 
derivative of the scalar differential group delay, the frequency derivative being a scalar 
second-order polarization mode dispersion function. The method further includes calculating 



a first result, according to 



d\f\ 




do 





, and multiplying a proportionality coefficient B2 by the 



first result. The mean differential group delay is calculated, therefore, in accordance with the 
following equation: 



(26) 
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[0017] Preferably, B2 is obtained from the probability density function of the scalar second- 
order polarization mode dispersion function. Most preferably, B2 is substantially equal to 
2.64. 



[0018] Yet another embodiment of the present invention provides a method for measuring a 
mean square differential group delay of at least one length of optical fiber, including an 
initial step of characterizing a polarization mode dispersion vector as a function of frequency 
using a frequency-domain polarization mode dispersion measurement apparatus. 
Additionally, a second-order polarization mode dispersion vector is calculated as a function 
of frequency from the polarization mode dispersion vector. The method further includes 
calculating a mean of the magnitude of the second-order polarization mode dispersion 
vector I to obtain a first result, according to ^|r^ | ^ . The first result is multiplied by a 

proportionality coefficient Aj to calculate the mean square differential group delay, in 
accordance with the following equation: 

^i(KI)=^Ls. (20) 

[0019] Preferably, Ai is obtained from the probability density fimction of the second-order 
polarization mode dispersion vector. Most preferably, Aj is substantially equal to 2.02. 

[0020] Still another embodiment of the present invention provides a method for measuring a 
mean square differential group delay of at least one length of optical fiber, including an 
initial step of measuring a magnitude of a polarization mode dispersion vector as a fimction 
of frequency, using a frequency-domain polarization mode dispersion measurement 
apparatus, where the magnitude of the polarization mode dispersion vector is a scalar 
differential group delay. The method further includes calculating a scalar second-order 
polarization mode dispersion function, the scalar second-order polarization mode dispersion 
function being a frequency-derivative of the scalar differential group delay. Additionally, a 

d\f\ 



first result is calculated according to 



dm 



The first result is multiplied by a 



proportionality coefficient, 5/, to calculate the mean square differential group delay, in 
accordance with the following equation: 
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[0021] Bi is preferably obtained from the probability density function of the scalar second- 
order polarization mode dispersion function. Bj is most preferably substantially equal to 
6.80. 

[0022] The measurement of polarization mode dispersion through an optical fiber may be 
through a single optical fiber link or a whole optical fiber route in an optical 
telecommunication network. 

[0023] In an additional embodiment of the method of the present invention, a method is 
provided for characterizing polarization mode dispersion of at least one length of optical fiber 
using a source of bandwidth B, The method includes an initial step of collecting polarization 
mode dispersion data as a function of frequency from a frequency-domain polarization mode 
dispersion measurement apparatus. The method further includes extracting one of a vector 
and a scalar frequency-dependent function from the polarization mode dispersion data, by 
applying a frequency-domain polarization mode dispersion technique, the one of the vector 
and the scalar function being one of a first-order and second-order polarization mode 
dispersion function. A systematic correction is applied thereto, wherein the systematic 
correction minimizes a systematic error in conventional methods of measuring mean 
differential group delay, caused by the finite source bandwidth B. Applying the systematic 
correction results in a derivation of one of a mean differential group delay (r) and a mean 
square DGD . 

[0024] Another method of measuring a mean differential group delay (r) of a length of 

optical fiber according to the present invention includes the steps of deriving a first mean (r) 

in accordance with equation (21), and a second mean (r) in accordance with equation (26), 

and then deriving a linear equation of a sum of the first mean (r) and the second mean (r) to 

calculate a combined mean (r) . A sum of coefficients of the linear equation is substantially 
equal to one. 
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[0025] Yet another method of measuring a mean square differential group delay of a 
length of optical fiber in accordance with the present invention includes the steps of deriving 
a first mean square differential group delay in accordance with equation (20), deriving 

a second mean square differential group delay in accordance with equation (25), and 
deriving a linear equation of a sum of the first mean square differential group delay r^j^ and 
the second mean square differential group delay to calculate a combined mean square 

differential group delay . A simi of coefficients of the linear equation is substantially 
equal to one. 

[0026] Other objects and features of the present invention will become apparent from the 
following detailed description considered in conjxmction with the accompanjdng drawings. It 
is to be understood, however, that the drawings are designed as an illustration only and not as 
a definition of the limits of the invention. 



Brief Description of the Drawings 

[0027] FIG. 1 is a pictorial view of an optical pulse propagating a length of birefringent 
optical fiber. 

[0028] FIG. 2 is a flow diagram illustrating a prior art method for calculating mean 
differential group delay (DGD). 

[0029] FIG. 3 is a block diagram illustrating a conventional interferometric polarization 
mode dispersion (PMD) measuring device. 

[0030] FIG. 4 is a block diagram illustrating a conventional apparatus used in a firequepcy- 
domain PMD measuring device. 

[0031] FIG. 5 is a plot of the standard deviation of the distribution of simulated mean DGD 
measurements of (r)^ , (|f^ |^ ) , and 
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[0032] FIG. 6 is a plot showing the improved standard deviation obtained with a systematic 
correction applied according to the method of the present invention. 

[0033] FIG. 7 is a plot of the standard deviation of the distribution of simulated square mean 
DGD measurements of (r ^ )^ , ^, (|f^ , and i(| 

[0034] FIG. 8 is a plot of standard deviation resulting from combining the methods shown in 
FIG. 5. 

Detailed Description 

[0035] The present invention relates to improving the accuracy of calculating a true mean 
differential group delay (DGD) or a true mean-square DGD (r) from a measurement of 

polarization mode dispersion (PMD). The present invention provides a method for correcting 
a systematic error in a conventional conversion of measured mean-square DGD 
measurements to a true mean DGD (r) . The systematic error stems from the use of a 
measurement source of a finite bandwidth B, and from the conventional conversion applied to 
measured root mean square DGD ^(r^)^ to obtain an estimate of tme mean DGD (r) . By 
incorporating improved averaging techniques to the measured root mean square DGD 
.J(r^ obtained with conventional PMD measurement apparatus, the present invention 
offers improved resolution in true mean (r) or true mean-square DGD <r^> measurements. 

[0036] The present invention also provides a method to improve accuracy in measurements 
of the true mean DGD (r) of optical fiber using statistics of a second-order polarization 
mode dispersion (SOPMD) and incorporating improved averaging techniques. A 30% 
reduction in measurement imcertainty is obtained when using the method in accordance with 
the present invention, which is equivalent to doubling the bandwidth B of an optical source, 

[0037] The following notation and definitions are used herein. Optical frequency (o refers to 
Incfk, where X is the wavelength of the source and c is the speed of light. Optical wavelength 
and frequency are used interchangeably to characterize the measurement source. 
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Analogously, spectral bandwidth and frequency bandwidth both refer to the optical source 
characteristics in wavelength and frequency space, respectively. 



[0038] An average over a finite angular frequency bandwidth 5 of a variable X(q)) associated 
with a given fiber is denoted and defined as follows: 



A "measurement of the mean" refers to a finite bandwidth average, as described by equation 
(1), applied to an experimentally measurable quantity. 

Applying equation (1), a root mean square of a parameter X is then defined as: 



A root mean square (RMS) of the DGD measured as a function of frequency over a source 
of and width B is therefore denoted generally by: 



where subscript "5 " indicates that measurements are taken over a finite source bandwidth, B, 

On the other hand a "true" RMS DGD, obtained by measuring over infinite bandwidth in 
equation (2) is denoted herein as . 

[0039] Equation (3) may be written as an integral of a function t^{(o) in accordance with 
equation (2). As well understood by those skilled in the art, however, when measuring 
polarization mode dispersion in the frequency domain, a DGD (r) is measured as the time 

delay between the principle states of polarization (PSP) associated with a discrete frequency, 
or equivalently, as the magnitude of the PMD measured at each discrete frequency over the 
source bandwidth 5. An RMS DGD value over a finite bandwidth B may therefore be 
calculated as an average of the discrete frequency-dependent DGD values over the source 
bandwidth, in accordance with the following: 



(1) 





(3), 
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where A<w = B/(2N+1X and a total of 2A^ + 1 measurements are made over source bandwidth 
B. 



[0040] A mean over an ensemble of N fibers is denoted by a bar over the parameter averaged, 
such as Xin the following equation: 

^4t^.. (5). 

Averaging of the same frequency-dependent parameter, the DGD for example, over an 
infinite number of fibers (N oo) with identically distributed birefiingence, where the 
parameter is measured for each fiber over a finite bandwidth 5, is equivalent to measuring the 
same frequency-dependent parameter of one fiber over an infinite bandwidth (B oo). For 
this "real" or "true" mean, taken over the limit of infinite bandwidth, the subscript is 
omitted. 

[0041] A second-order polarization mode dispersion (SOPMD), which is the derivative of the 
polarization mode dispersion vector with respect to frequency, is represented by . 

[0042] A fiber link may refer to any fiber optic cable or cable assembly connecting two 
points, such as optical amplifiers, along an optical transmission line. A fiber route may 
include multiple fiber links, and includes the total fiber optic cable between any two 
terminating points in an optical transmission line. 

[0043] In a first embodiment, a method is provided to more accurately measure the mean of 
the differential group delay from conventional PMD measurement techniques. The method 
employs a systematic method to minimize a systematic error that occurs when converting 
root-mean-square measurements of the spectral-averaged DGD to the true mean DGD (r) 

using conventional means. The error is a consequence of the finite bandwidth B of the 
measurement source. 

[0044] Referring to FIG. 2, a conventional method 20 of estimating the true mean DGD (r) 
includes the step 22 of characterizing the PMD vector using techniques commonly known, 
such as Jones Matrix Eigenanalysis (JME), MuUer Matrix Method (MMM), and Poincare 
Sphere Analysis (PSA). In step 23, the mean DGD is calculated as an average over the 
spectral source bandvsddth to estimate a true mean. Alternatively, in step 24, the RMS value 
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of the DGD i^J{^^)^ is calculated from the magnitude of the frequency-averaged PMD 
vector, in accordance with equation (4), for example. Altemately, step 26 is used to derive 
^(t^ directly from an interferometric measurement of the square of time delay averaged 

over the source bandwidth. In step 28, ^J{^^ derived from either step 24 or step 26 is 
multiplied by a conventional correction factor to calculate an estimate of tme mean (r) 29. 

[0045] It has been conmionly assumed in the prior art that the DGD can be statistically 
represented by a Maxwellian distribution. Therefore, based on statistical probability theory, 

an estimate of mean DGD is calculated by measuring ^^{t^)^ for a fiber and estimating the 

mean that would be found for an ensemble comprised of an infinite niunber of fibers in 
accordance with equation (5), by multiplying the result by the conventional correction factor 

of /— . In other words, by conventional methods, the following equation is applied to 
V 3^ 

obtain an estimate of true mean (r) 29 in accordance with step 28 of FIG. 2: 



(6) 

The validity of this approximation may be verified by comparing: 



where the ensemble averaging is taken over a multiplicity of different fibers with identically 
distributed birefiingence. 

[0046] Equation (7) is only correct, however, as 5 ^ oo , and introduces a systematic error 
when the measurement bandwidth B is finite. This is the easiest to see in the extreme case of 
a very small bandwidth, that is as 5 -> 0 . This case corresponds to a situation where only one 
frequency point per fiber is used to make a measurement of mean DGD. Averaged over 
many fibers with identically distributed birefiingence, the left hand side LHS of equation (7) 
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will converge to (r) , while its right hand side RHS will converge to /— (r) , resulting in 

V 3;r 

about an 8% systematic error. 



[0047] In the method of the present invention, a function F is derived, which is defined on 
the frequency span from -8/2 to B/2 to obtain an estimate of the DGD as a function of source 
bandwidth 5, denoted as '^{b) , such that 

F\t{co)-\=t{b), (8) 

The function is defined so that for an ensemble of many different fibers with the same , 
its ensemble average converges to the true mean (z") , as 



(9) 



and the standard deviation of r{B) , denoted as stdijiB)) , is minimized relative to the 
standard averaging. 

[0048] To derive a systematic correction factor, the systematic error is estimated 
asymptotically when the product of the bandwidth B and is much greater than 1 

( ^RMs^ » 1 ) and the systematic error is small. The result of any measurement of the mean 
square DGD averaged over a finite source bandwidth 5, or ^ , is represented by the 
foUowdng equation: 

^(^')^=W'+'7. (10) 

where r] represents the deviation of the measurement from the mean found for the 
ensemble, and the average of ri values for the ensemble is zero (7=0). Assuming that this 
deviation is small, a binomial expansion can be applied to approximate 7 so that equation 
(10) can be written in the following form: 
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[0049] Since 7=0, taking an ensemble average of each side of equation (11) leaves only the 
rf^ term and (r) on the right hand side, so that the systematic error s in measuring the mean 
DGD can be estimated from the following equation: 



= (r)-l±/r^\ =jf—. (12) 



8(r) 

[0050] The systematic error s is equivalently a systematic correction factor s that is applied 
to the measured RMS DGD ^(r^)^ to improve the resolution of true mean DGD (r) 
calculations. 

[0051] The variance of {r^) ^ was evaluated in M. Shtaif and A. Mecozzi, "Study of the 

Frequency Autocorrelation of the Differential Group Delay in Fibers with Polarization Mode 
Dispersion," IEEE Photon, Teck Lett, Vol. 25, pp. 707-709 (May 2000) for r^B » 1 as 
follows: 

var(r^)^=(;rVV2)(r)VB, (13) 
where the variance is defined as follows: 

H^')A^').-{^')1 '''' 

[0052] With the help of equation (10) to derive an expression for 7^ , and substituting the 
expression for the variance from equation (13), a value for the systematic error can be 
evaluated from equation (12) as follows: 
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[0053] In order to obtain an unbiased estimate of the mean DGD (r) , therefore, the 
following equation is preferably applied to the mean-square DGD (r^)^ measured using any 
of the conventional measurement techniques for obtaining the mean-square DGD: 



"""" ' ° ■ 

[0054] One embodiment of the method of the present invention, therefore, includes the steps 
of deriving a frequency-averaged RMS measurement of the DGD i^^)^ from conventional 
PMD measuring apparatus, and applying a systematic correction to the measured RMS to 
obtain true mean DGD (r) . As shown in equation (16a), the systematic correction is made 
by adding a systematic correction factor e to the conventional approximation of (r) that is 
described by equation (6). The preferred systematic correction factor e is characterized by 
equation (16b) above. The measviring apparatus can include the time-domain and frequency- 
domain measurement devices, as described below. 

[0055] Examples of conventional measurement techniques to which this systematic 
correction can be applied include, but are not limited to, frequency-domain measurement 
techniques such as JME, MME, PSA, and time-domain techniques such as low-coherence 
interferometry. A fimdamentally equivalent method to a scaiming Michelson interferometer, 
that also works in the frequency domain, is the Fixed Analyzer Technique. 

[0056] A block diagram of an interferometer 30 for PMD measurements in the time-domain 
is shown in FIG. 3. Time-domain techniques provide a direct measure in the time-domain of 
the frequency-averaged DGD. The commonly known interferometric technique for 
measuring PMD is described in more detail in N. Gisin, J.P. Von der Weid and J.P. Pellaux, 
"Polarization Mode Dispersion of Short and Long Single-Mode Fibers," J. Lightw. TechnoL, 
Vol. 9, p. 821 (1991), which is incorporated herein by reference. In the Michelson-type 
interferometer 30 of FIG. 3, for example, a spectrally broad (low-coherence) source 32 is 
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used to couple light 34 into the device under test (fiber) 36. The output of fiber 36 is 
coUimated using coUimating optics 38 and directed onto a beamsplitter 40. The beamsplitter 
40 divides the light into two paths along two aims of the interferometer 30. One path is 
directed toward a moveable mirror 42, which can be translated along a path coincident with 
the direction of the incident light. The second is directed toward a stationary mirror 44. 
When the mirror 42 is translated, interference fringes are seen by a detector 46 as long as the 
time-delay difference between the two arms matches a delay generated in the device under 
test 36 to within a coherence time of the source 32. 

[0057] In the simplest case of a non-mode-coupled device, the interferogram provided by a 
conventional interferometer 30 is symmetric and represented by two side lobes and a central 
autocorrelation peak. The distance from a side lobe to the central peak represents the DGD 
averaged over the source spectrum or bandwidth. In the mode-coupled case of a fiber, the 
interferogram includes muhiple peaks representing various principle states of polarization 
(PSP) associated with different frequencies. In this case, as well-known to those skilled in 

the art, the DGD is characterized as a RMS value over the bandwidth, that is ^{^^) ^ , as 

described by equation (3). Discrete measurements of the PMD as a function of frequency are 
not needed since the time-domain measurements are inherently frequency-averaged over the 
sovu-ce bandwidth. 

[0058] FIG. 4 shows a block diagram of a common polarimeter 50 that can be used in a 
measurement device for any of the frequency-based techniques, such as PSA, JME, and 
MMM. In FIG. 4, a polarization state of a tunable laser 52 is controlled by a polarization 
controller 54. The polarized output from the controller 54 is coupled to the fiber under test 
56. The polarization state of the optical signal after traversing the length of fiber 56 is then 
analyzed with a conventional polarization state analyzer 58. By stepping the laser 52 in 
specified frequency intervals over the entire bandwidth, an output polarization state for each 
input polarization state is obtained as a fimction of frequency. Each technique differs in 
which input polarization states are used and in how the data set of measured output 
polarization states are manipulated. The result, however, is the same: a measurement of the 
PMD vector, that is, both DGD and the PSP as a fimction of frequency. 

[0059] The PSA and JME techniques are well-known to those skilled in the art and are 
described in detail, for example, in Normand Cyr, Andre Girard, and Gregory W. Schinn, 
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"Stokes Parameter Analysis Method, the Consolidated Test Method for PMD 
Measurements", NFOEC '99 Convention, Chicago, III (1999), which is incorporated herein 
by reference. These techniques are based on the measurement of so-called polarimetric 
Stokes parameters. These two techniques differ in the algorithms used to describe the PMD 
from the measured Stokes parameters, as well as in the required form of the generally three 
measured Stokes parameters. The JME algorithm requires only linear input states of 
polarization (SOP), whereas the PSA algorithm usually requires one of the three input states 
to be circular. 

[0060] The PSA and JME methods are alike, however, in requiring the same test procedures. 
Both find the PSP and the DGD as a fimction of frequency from raw polarimetric data. The 
raw polarimetric data is the normalized Stokes vector representing the SOP at the output of a 
fiber as a fimction of optical frequency for a number of input SOPs. The normalized Stokes 
vector, s, is the position of the output SOP on a Poincare sphere. Basically, both methods 
measure the local "angular velocity" (where time is actually optical frequency here) of the 
output SOP from a measurement of its postion on the sphere at two neighboring optical 
frequencies o) and (co+Aco). One difference is that PSA is performed directly in the Stokes 
parameter representation of the SOP, while JME is performed in the Jones vector 
representation, after the Stokes vectors have been transformed into normalized Jones vectors. 

[0061] The MMM technique has been described as a combination of the PSA and JME 
method. The MMM method is described in detail in R.M. Jopson, L.E. Nelson, and H. 
Kogelnik, , "Measurement of Second-Order Polarization-Mode Dispersion Vectors in Optical 
Fibers", IEEE Phot, Tech, Lett., Vol. 1 1, No. 9, pp. 1 153-55 (Sept. 1999). Like the other two 
techniques, the MMM requires the determination of transmitted output Stokes vectors 
corresponding to two input polarizations s, at each frequency to. Unlike PSA however, 
MMM has the ability to determine the rotation matrix of the fiber. It also stays entirely in 
Stokes vector space. The technique is advantageously independent of the relative angle 
between the two input linear polarizations. MMM may also be used with a interleaving 
method to measure the first- and second-order PMD vectors of fibers, as described in Jopson, 
etal. 

[0062] Both the interferometer of FIG. 3 and the polarimeter of FIG. 4, as well as any 
conventional PMD apparatus, can be used to produce frequency-averaged RMS 
measurements of the DGD over the source bandwidth. The systematic correction obtained 
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above in equation (16b) can then be directly applied to the measured RMS DGD Tj^^ as in 
equation (16a) to minimize the error in calculating the true mean DGD (r) . 



[0063] The method of the present invention provides a systematic correction for the property 



that measurements of (r)^ and n/— (^^ ) are not equivalent, due to the physical limitation 



imposed by the finite source bandwidth. In other words, taking the square root of a 
frequency-averaged parameter leads to systematic errors for finite bandwidth measurements. 
The systematic correction for this property can also be extended to higher powers of DGD 
and to frequency derivatives of PMD. 

[0064] In another embodiment of the present invention, a method of more accurately 
measuring mean DGD <r> is provided by measuring a second-order PMD (SOPMD) over a 
finite bandwidth B and applying a systematic correction to minimize the random errors that 
typically occur from the measurement over a finite bandwidth. The systematic correction 
includes an averaging procedure that avoids the systematic error induced by conventional 
methods of taking roots after averaging over fi-equency. This method produces superior 
results when the fiber length over which measurements are performed is considerably longer 
than the birefringence correlation length, as is the case for fiber spans in teleconununication 
systems. 

[0065] It is shown in G.J. Foschini and CD. Poole, "Statistical Theory of Polarization 
Dispersion in Single Mode Fibers," lEEEJourn, of Lightwave Techn., Vol. 9, pp. 1439-1456 
(Nov. 1991) ("Foschini and Poole"), that the mean of the SOPMD in a fiber, represented by 
{^a, ) depends only on its mean square DGD (t^ ) when the fiber is much longer than the 
correlation length of its birefiingence, in accordance v^th the following equation: 



[0066] It has also been shown that the SOPMD can be accurately measured experimentally if 
the measurements are performed with a small enough frequency step. {See R. M. Jopson, L. 
E. Nelson, H. Kogehiik, "Measurement of Second-Order Polarization-Mode Dispersion 
Vectors in Optical Fibers", IEEE Phot Tech, Lett, Vol. 11, pp. 1153-55 (Sept. 1999), 
incorporated herein by reference. ) Jopson et al. applied an interleaving technique to increase 
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the frequency resolution of the measurements, then applied the MuUer Matrix Method 
(MMM) to these measurements to obtain the PMD vector, and calculated the RMS DGD 
^RMs from the magnitude of the PMD vector f . In addition, Jopson et al. calculated the 
SOPMD from these measurements and applied the Foschini and Poole relation described by 
equation (20) to estimate either the mean square DGD ^r^^ or the root mean square DGD 

[0067] As described above, however, due to the finite bandwidth B of the source used for the 
measurement, a systematic error occurs when taking the square root of a mean square value 
averaged over frequency. The direct application of equation (17) according to the prior art, 
therefore, produces inaccurate results, especially over low source bandwidths. The error is 
present whether it is a first, second, or higher order of PMD that is measured as a frequency- 
averaged RMS value. The estimate of the mean DGD (r) from second-order PMD 

measurements, therefore, will also be biased at small measurement bandwidths for reasons 
very similar to the embodiment described above with respect to first order PMD 
measurements, 

[0068] An additional embodiment of the method of the present invention, therefore, provides 
an accurate measure of the mean DGD (r) from measurements of the SOPMD taken over a 
finite bandwidth jB, by applying a systematic correction to minimize the random errors that 
occur from the measurement over a finite bandwidth. The systematic correction includes an 
averaging technique that avoids the systematic error induced by conventional methods of 
taking roots after averaging over frequency. 

[0069] A further embodiment of the method of the present invention provides an accurate 
measure of the mean squared DGD from measurements of the SOPMD taken over a 
finite bandwidth, by applying a systematic correction to minimize the random errors that 
occur from the measurement over a finite bandwidth. The systematic correction includes an 
averaging technique that avoids the systematic error induced by conventional methods of 
taking roots after averaging over frequency. 

[0070] The SOPMD is the first derivative of the PMD with respect to frequency, representing 
the change in the PMD as a function of frequency. Specifically, if rico) denotes the PMD 
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vector for a frequency q>, defined by the principal states of polarization (PSP) at that 



frequency co, then the magnitude of the vector 



is the differential group delay DGD r 



or the set of PSP defined for frequency co. Therefore, a set of vector PMD data as a function 
of frequency © and a set of scalar DGD data as a function of frequency © may be generated 
from any of the methods used to extract the PMD from a frequency-domain measurement 
device. Such measurement device may comprise a polarimeter or a Fixed Analyzer, for 
example. 

[0071] In addition, a second-order PMD can be defined as either a vector or scalar quantity. 
The vector SOPMD is denoted as and is found from the set of vector PMD data, 
according to the conventional method of defining a derivative from a discrete data set of n 
sampled frequency points as follows: 



(18) 



[0072] The scalar SOPMD is more easily calculated as the derivative of the scalar DGD and 



d 

is denoted as 



, where again, the magnitude of the vector 



is the differential 



dco 

group delay DGD r , a scalar quantity. This scalar SOPMD is found directly from the set of 
frequency-dependent scalar DGD t values, i.e. from r(<u) , according to the conventional 
method of defining a derivative from a discrete data set as follows: 



(19) 



[0073] Different quantities scaling as r and respectively are preferably used to avoid bias 
in the calculated results for each of these embodiments. A probability density ftinction (PDF) 
for the magnitude of the scalar SOPMD and the vector SOPMD have been calculated by 
Foschini, et al. in "Probability Densities of Second-Order Polarization Mode Dispersion 
Including Polarization Dependent Chromatic Fiber Dispersion," IEEE Phot. Tech. Lett., Vol. 
12, pp. 293 -295 (March 2000), (Foschini et al.), which is incorporated herein by reference. 
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[0074] By integrating the probability density function for the magnitude of the vector 
SOPMD in conjunction with equation (17), the mean square and the mean of the DGD may 
be derived, resulting in the foUovsdng relationships: 

A{K\) = ^iMs (20) 

9 

= (21) 

«1.37 (Catalan constant 

G«0.916) may be calculated from the PDF of the SOPMD, and verified by simulations. 
Written this way, equations (20) and (21) provide a way to estimate, without bias, the mean 
DGD and the mean square DGD using the SOPMD vector magnitude \f^\ . 

[0075] PMD and SOPMD measurements are strongly correlated as seen from the relation for 
the conditional mean square of the SOPMD, derived in J.P. Gordon, "Statistical Properties of 
Polarization Mode Dispersion", Proc. VSSV2 (Venice, Italy 2002) and provided below: 



Coefficients A, — 



^ » 2.02 and 



(-1)' 



{^■iH^vHi^'Y (22) 

Therefore, if a specific measurement of the mean DGD of a given fiber is above the "true" 
mean, the mean DGD estimated from equation (21) will be overestunated as well. Gordon, 
however, has shown that the conditional mean square of the component of the second order 
PMD vector that is parallel to the direction of the PMD vector is independent of the PMD 
vector length and can also be expressed in terms of the mean DGD (r) as follows: 

(23) 

[0076] Equation (23) can also be used to provide an independent measiu-e of the mean square 
DGD (r ^ ) by recognizing that the parallel component of the second order PMD in equation 
(23) is, equivalently, the frequency derivative of the DGD, as follows: 

r.-^ • (24) 
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[0077] The frequency derivative may be calculated from the scalar DGD vs. frequency 
measurements performed by a standard non- vectorial instrument, allowing the technique to 
be used with conventional PMD test sets. As before, to avoid systematic errors, equation (23) 
may be combined with equation (24) and rewritten after application of the PDF of the scalar 
SOPMD as: 



din 



dCD 



(25) 



and 



^-1 



(26) 



These relationships provide another imbiased evaluation of the mean DGD using the 
frequency derivative of the DGD, and can be verified with simulations. The coefficients 

-1 



5,= 



« 6.80 and 5, = 
21og2 ' 



1 



« 2.64 do not depend on the fiber 



parameters, and may be calculated from the PDF derived for SOPMD, as described in 
Foschini et al. 

[0078] While equations (20) and (21) require the SOPMD in vector form, found in 
accordance with equation (18) above, equations (25) and (26) require only the scalar DGD, 
measured at discrete frequency intervals, so that the scalar SOPMD may be calculated, for 
example, in accordance with equation (19). Non- vectorial instruments that provide discrete 
rather than frequency-averaged scalar DGD measurements may, therefore, be used to 
calculate r^^^ and (r) from equations (25) and (26) respectively. 

[0079] The improved accuracy of the method in accordance with the present invention for 
estimating the mean DGD (r) and the mean square DGD (r^) from the measured SOPMD, 
to that of the conventional averaging techniques, has been demonstrated in simulations. 

[0080] To estimate the measurement accuracy numerically, a simulation may be performed 
on 10000 fibers consisting of 1000 random sections of birefringence. The polarization 
components, for example. Stokes components as well-known in the art, of the birefringence 
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vector of each section in the simulation preferably have independent identical Gaussian 
distribution. The delay of all birefringent sections is preferably normalized to ensure that the 
true mean DGD is equal to Trms, and to prevent additional variations in the measurement due 
to the difference in "true" mean DGD among fibers. The PMD vector and its frequency 
derivative are preferably computed at 256 angular frequency points separated by 
In 

A<a) = , which is much less than the DGD correlation bandwidth. This corresponds to 

the DGD-angular bandwidth product Bxj^j^ « 150 . By averaging DGD and the other 
quantities in equations (20), (21), (25), and (26) over a different number of frequency points 
for every fiber, the effect that the bandwidth of the optical source has on the accuracy of the 
PMD measurement can be evaluated. It is expected that all simulated measurements will be 
centered at (r) independently of the measurement bandwidth used. 

[0081] FIG, 5 shows a plot of the normalized standard deviations of the mean DGD 
measurements obtained by direct averaging of DGD (solid line) 70, by estimating (r) from 
the vector SOPMD using equation (21) (dashed line) 72, and from the scalar SOPMD using 
equation (26) (dotted line) 74. As shown in FIG. 5, the mean DGD estimated from the 
second order PMD, whether estimated from the vector or scalar SOPMD, provides 
significantly better accuracy over the prior art method of direct averaging. 

[0082] Referring to FIG. 6, reductions in the standard deviation achieved by applying the 
statistical SOPMD to estimates of the mean DGD are shown. In FIG. 6, the number of 
samples 77 (pps) is plotted as a fimction of the mean DGD calculated from the vector 
SOPMD 75, and as a fimction of the mean DGD calculated directly as an average 79. The 
standard deviation for the averaged DGD 79 is about .044 ps, as compared with a standard 
deviation of .028 pps for tke mean DGD derived from SOPMD. This amounts to a factor of 
about a 42 reduction in the standard deviation, which corresponds to a two-fold reduction of 
the measurement variance. This reduction is comparable to the effect that would be achieved 
by doubling the normalized measurement bandv/idth Bxj^ in standard PMD measurement 
techniques. 

[0083] FIG. 7 shows a plot of the normalized standard deviations for the mean square DGD 
^^)b ' directiy (solid line) 80, and by estimating the mean square DGD from the 

vector SOPMD using equation (20) (dashed line) 82, and from the scalar SOPMD using 
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equation (25) (dotted line) 84. A significant reduction of the measurement uncertainty of at 
least 30% is shown in the simulation plotted in FIG. 7. 

[0084] These numerical simulations also show that the mean DGD measurements taken by 
different methods are not entirely correlated. Averaging the results of these estimation 
techniques, therefore, allows further reduction of the measurement uncertainty. FIG. 8 
presents the incremental improvements in the measurement accuracy obtained by combining 
measurement techniques. Shown in dashed cxirve 90 is the standard deviation of the 

measurements of y^2^|f^|^^ = (r) , which represents the preferred method of calculating (r) of 

those shown in FIG. 5. The dotted curve 92 represents the standard deviation resulting from 
combining the estimates of DGD (r) from the vector SOPMD 72 and the scalar SOPMD 74 
(see FIG. 5), Marginal improvement in the standard deviation is shown by combining 
techniques as shown by dotted curve 92. The solid curve 94 additionally averages in the 
result from direct averaging of the DGD 70 (see FIG. 5), which also shows some 
improvement over the use of the single technique 90. 

[0085] An alternate method of the present invention includes combimng the statistical results 
of more than one of the methods provided herein to further reduce the error in measuring the 
DGD. Specifically, a linear combination for the mean differential group delay DGD (t) is 
preferably obtained from the LHS of equations (21) and (26), to obtain a measurement of 
mean differential group delay DGD (r) . Coefficients of each term in the linear combination 
are preferably normalized so that a sum of the coefficients is substantially equal to one. 
Similarly, a linear combination for the mean square differential group delay DGD ^ ) is 

preferably obtained from the LHS of equations (20) and (25), to obtain a measurement of 
mean square differential group delay DGD ^r^)^ . Coefficients of each term in the linear 

combination are preferably normalized so that a sum of the coefficients is substantially equal 
to one. Analyses of the statistics of even higher orders of the PMD may also offer 
improvement over the prior art. 

[0086] The methods provided herein in accordance with the present invention may be applied 
to more accurately measure the true mean DGD in a fiber link or a whole fiber route in an 
optical telecommunication system. 
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[0087] Although illustrative embodiments of the present invention have been described 
herein with reference to the accompanying drawings, showing results of applying the 
statistics up to the second-order PMD to more accurately measure the differential group delay 
of fiber, it is to be understood that the invention is not limited to those precise embodiments, 
and that various other changes and modifications may be effected therein by one skilled in 
the art without departing firom the scope or spirit of the invention. 
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